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O ' We present deep mid-infrared spectroscopy with the Spitzer Space Telescope 

of 13 submillimeter galaxies (SMGs) in the Great Observatories Origins Deep 
Survey North (GOODS-N) field. We find strong polycyclic aromatic hydrocar- 
Q_i' bon (PAH) emission in all of our targets which allows us to measure mid-IR 

O ■ spectroscopic redshifts and place constraints on the contribution from star for- 

mation and active galactic nuclei (AGN) activity to the mid-IR emission. In 
the high signal-to-noise ratio composite spectrum of SMGs, we find that the 
hot dust continuum from an AGN contributes at most 30% of the mid-IR lu- 
minosity. Individually, only 2/13 SMGs have continuum emission dominating 
the mid-infrared luminosity; one of these SMGs, CI, remains undetected in the 
. deep X-ray images but shows a steeply rising continuum in the mid-IR indica- 

tive of a Compton-thick AGN. We find that the mid-IR properties of SMGs are 
distinct from those of 24/im-selected ULIRGs at z ~ 2; the former are predom- 
inantly dominated by star formation while the latter are a more heterogeneous 
sample with many showing significant AGN activity. We fit the IRS spectrum 
and the mid-IR to radio photometry of SMGs with template spectral energy dis- 
tributions to determine the best estimate of the total infrared luminosity from 
star formation. While many SMGs contain an AGN as evinced by their X-ray 
properties, our multi-wavelength analysis shows that the total infrared luminos- 
ity, Lxr, in SMGs, is dominated by star formation and not AGN activity. We 
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ffiid that high redshift SMGs lie on the relation between L IR and Lpah,6.2 (or 
-kpAH,7.7 or -^pah,h.3) that has been established for local starburst galaxies. This 
suggests that PAH luminosity can be used as a proxy for the star formation rate 
in SMGs. SMGs are consistent with being a short lived cool phase in a massive 
merger. Although there appears to be an AGN present in many SMGs, it does 
not appear to have become strong enough to heat the dust and dominate the 
mid- or far-infrared emission. 

Subject headings: galaxies: evolution — galaxies: starburst — galaxies: active 
— infrared: galaxies — submillimeter — techniques: spectroscopic 

1. Introduction 

The COsmic Background Explorer (COBE) revealed that the extragalactic background 
light at infrared (IR) wavelengths was comparable to that at optical wavelengths indicating 
that half of the star formation and active galactic nuclei (AGN) emission, averaged over 
cosmic time, is obscured by dust (e.g. Puget et al. 1996; Hauser & Dwek 2001). In the local 
Universe, the volume density of dusty, infrared luminous galaxies is insufficient to account for 
this background. Therefore the density of such galaxies must evolve strongly with redshift 
(e.g. see Lagache, Puget, & Dole 2005 for a review). In fact, imaging observations with the 
Infrared Space Observatory (ISO) and the Spitzer Space Telescope (hereafter Spitzer) have 
shown that about 70% of the comoving star formation rate density (SFRD) at 0.5 < z < 3 
is obscured by dust (Chary & Elbaz 2001; LeFloc'h et al. 2005; Perez-Gonzalez et al. 2005). 

In the local Universe, we refer to these dusty galaxies as luminous infrared galaxies 
(LIRGs, 11 < log (Lir [L Q ]) < 12) and ultraluminous infrared galaxies (ULIRGs, log (Lir [L ]) > 
12) where Lir is the total IR luminosity from 8-1000 /iin (see Sanders & Mirabel 1996 for a 
review). It is now becoming clear that IR luminous galaxies are a heterogenous population 
with selection functions at different mid- and far-infrared wavelengths revealing different 
subsets of this population (Daddi et al. 2005; Lutz et al. 2005b; Yan et al. 2005; Pope et 
al. 2006). 

Perhaps the most enigmatic members of the high redshift ULIRG population are the 
submillimeter (submm) galaxies (SMGs). These galaxies were initially discovered in the late 
1990s with the Submillimetre Common User Bolometer Array (SCUBA, Holland et al. 1999) 
on the James Clerk Maxwell Telescope (JCMT) and then later with MAx-Planck Millimetre 
BOlometer (MAMBO) camera on the Institut de Radio Astronomie Millimetrique (IRAM) 
30-m telescope (see Blain et al. 2002 for a review). They have proven to be difficult to study, 
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due to the large beam size of submm observations and the difficulty in finding counterparts 
at other wavelengths (e.g. Ivison et al. 2002; Pope et al. 2005). Nevertheless, it has been 
found that they are very massive systems (Swinbank et al. 2004; Borys et al. 2005; Greve 
et al. 2005; Tacconi et al. 2006) at z ~ 2 (Chapman et al. 2005; Pope et al. 2006; Aretxaga 
et al. 2007) which are thought to be connected to the most massive present day galaxies via 
an evolutionary sequence (Lilly et al. 1999). A crucial aspect that needs to be explored in 
SMGs is what powers their extreme luminosities: AGN activity or star formation? 

The mid-IR spectral regime is a particularly good probe of what is producing the total 
infrared luminosity since it has contributions from both star-formation and AGN each with 
a unique signature at these wavelengths. The 5-15 /zm mid-IR spectral energy distribution 
(SED) can be decomposed into 3 main components (e.g. Sajina et al. 2007): (1) emission 
features from what are thought to be polycyclic aromatic hydrocarbon (PAH) molecules 
(Puget & Leger 1989; Allamandola et al. 1999); (2) power-law (or warm blackbody) emis- 
sion; and (3) extinction characterized by prominent silicate absorption features (e.g. Draine 
2003). The first is understood to be powered entirely by star formation while the second 
is predominantly emission from the AGN (Genzel et al. 1998). The second component can 
also contain contributions from hot dust (very small grains) present in the most energetic 
HII regions (e.g. Tran et al. 2001). The third component can be found in both starbursts 
and AGN since it just requires a dust screen around hot dust emission. 

Measuring the contribution from ULIRGs to the global star formation rate requires dis- 
secting their IR emission into contributions from starbursts and AGN. Mid-IR spectroscopy 
with ISO made it possible to decompose the various contributors in local ULIRGs with 
the result that they are predominantly starburst powered, although they also show signs 
of AGN activity (e.g. Genzel et al. 1998; Lutz et al. 1999; Rigopoulou et al. 1999; Sturm 
et al. 2000; Tran et al. 2001). In particular, ISO studies of local ULIRGs showed that at 
low luminosities, the ULIRGs are powered mainly by a starburst but at the very highest 
luminosities (Lir ^3x 1O 12 L ), the dominant emission source is the AGN (Lutz et al. 1998; 
Tran et al. 2001). From this one might expect that SMGs, with their extreme IR luminosities 
(L IR ~6 x 1O 12 L , Pope et al. 2006), are primarily AGN-dominated. 

The Spitzer Infrared Spectrograph (IRS; Houck et al. 2004), with its sensitivity and 
spectral resolution, can extend this effort outside the local Universe, allowing much larger 
samples to be observed to much fainter flux limits. The IRS can detect emission features and 
continuum in ULIRGs out to z ~ 4 (Valiante et al. 2007). The 3 mid-IR spectral components 
listed above have been fit to the Spitzer IRS spectra of ULIRGs at various redshifts to 
determine the AGN contribution to the mid-IR luminosity (e.g. Sajina et al. 2007). Although 
X-ray imaging is very effective at identifying the presence of an AGN in a galaxy, it is 
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subject to obscuration in extreme, Compton-thick sources. The X-ray luminosity provides 
a direct measurement of the power of the AGN at X-ray energies but it is necessary to 
make assumptions about the geometry of the obscuring dust to determine the bolometric 
luminosity of the AGN (e.g. Alexander et al. 2005). Furthermore, AGN indicators at UV and 
optical wavelengths can be very obscured in dusty SMGs (Swinbank et al. 2004; Chapman 
et al. 2005). Mid-IR spectroscopy, by directly detecting the thermally re-radiated emission 
from the obscuring dust, provides one of the most effective ways to identify AGN, including 
Compton-thick sources, and measure their contribution to the infrared luminosity (e.g. Lutz 
et al. 2004). 

Two previous studies have presented IRS spectroscopy of SMGs. Menendez-Delmestre 
et al. (2007) presented IRS spectra of 5 radio-detected, spectroscopically confirmed SMGs. 
They found that four SMGs with lower redshifts (z < 1.5) showed PAH emission and a 
composite spectrum of these four sources fits well to a scaled M82 spectrum plus a power-law 
component. They conclude that these systems host both star formation and AGN activity. 
Their fifth source at z — 2.4 was found to have a more AGN-type spectrum. Valiante et 
al. (2007) discussed IRS spectra of 13 bright SMGs, 9 of which are detected with the IRS 
(the first two spectra from this sample were presented in Lutz et al. 2005a). This sample 
was more representative of the redshift distribution of SMGs. They also found that SMGs 
are likely to be mainly starburst powered. Sajina et al. (2007) used IRS spectra to study 
a sample of Spitzer 24 /zm-selected galaxies at z ~ 2. This sample was selected to have 
SWjum > 900 /iJy and very red 24-to-8/im and 24-to-i? colors (Yan et al. 2007). 75% of this 
sample consists of continuum-dominated sources, with some weak PAH emission indicative 
of star formation activity. They conclude that these sources display mid-IR spectra similar 
to those of the most luminous local ULIRGs. Follow-up studies of these galaxies have shown 
that they are not generally (sub)mm detected sources (Si. 2m m ^ 1.5 mJy for 33/40 sources, 
Lutz et al. 2005b). 

In this paper we present Spitzer IRS spectra of a sample of SMGs in GOODS-N which 
have extensive multi-wavelength information from X-ray to radio wavelengths. This allows 
us to produce full spectral energy distributions and to calculate the bolometric luminosity. 
Our Spitzer IRS observations go substantially deeper (S^^m > 200 \iJy) than previous 
work. This enables a study of the typical blank-field submm galaxy population which have 
£850 Atm. > 2 mJy. The main goal of our program is to decompose the mid-IR spectra into star 
formation and AGN components, in order to infer how much of the total infrared luminosity 
is coming from each. From the mid-IR spectroscopy we are also able to obtain spectroscopic 
redshifts and measure the luminosities and equivalent widths of individual PAH lines to 
investigate various line diagnostics. We compare our measurements with line strengths in 
other star-forming galaxies to put constraints on the evolution of SMGs and the relationship 
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between star formation and black hole growth. 

This paper is formatted as follows. Our Spitzer observations are described in Section 2 
and the data analysis is explained in Section 3. Section 4 presents the main results of the IRS 
spectroscopy including redshifts, spectral decomposition and AGN classification. In Section 
5 we fit the full SEDs of these SMGs to determine their total IR luminosities. We explore 
correlations between the PAH and IR luminosities in Section 6. We discuss our results in the 
context of the role that SMGs play in galaxy evolution in Section 7. Finally, our conclusions 
are listed in Section 8. We include an appendix which discusses the noise properties of our 
Spitzer IRS spectra and includes notes on individual sources. 

All magnitudes in this paper use the AB system. We use a standard cosmology with 
H = 73kms- 1 Mpc~ 1 , Q M = 0.3 and Q A = 0.7. 

2. Observations 

The Spitzer IRS observations presented here were taken as part of Spitzer program GO- 
20456 (PI: R. Chary). The full sample consists of high redshift SMGs, AGN and optically 
faint 24 /im sources in the Great Observatories Origins Deep Survey (GOODS, Dickinson et 
al. 2003) North field. The GOODS-N field is one of the most extensively studied regions 
of the sky, with deep data existing across all wavebands including: Chandra 2 Msec X-ray 
observations (Alexander et al. 2003); deep HST optical imaging using the Advanced Camera 
for Surveys (ACS) in four bands, -B435, Vqoq, ^775 and zg$o (Giavalisco et al. 2004); and deep 
Spitzer imaging at five infrared wavelengths (Dickinson et al., in preparation). More recently, 
GOODS-N has been surveyed with IRS peak-up observations at 16 /itn (Teplitz et al. 2005) 
and MIPS imaging at 70 /im (Frayer et al. 2006). In addition to the extensive space-based 
imaging campaign, many ground-based programs have also targeted this field for imaging 
and spectroscopy making it one of the best field for investigating the role of AGN and star 
formation in high redshift galaxies. Here we present the IRS spectroscopy for the SMGs 
from this project. Results from the other sources will be presented in other papers. 

Our sample of SMGs was primarily selected from sources detected in the GOODS-N 
SCUBA 'super-map' with secure mult i- wavelength counterparts (Borys et al. 2003; Pope et 
al. 2006; Pope 2007). We imposed a 24 /im flux cut of S24 > 200 /iJy in order to ensure that 
we achieve adequate signal to noise (SNR) on the continuum in the IRS spectra. Roughly 
half of the secure counterparts in Pope et al. (2006) meet this 24 /xm flux cut which yields a 
sample of 10 IRS targets (after removing sources which are already observed with the IRS in 
other programs). To expand our sample slightly, we also included three SCUBA sources from 
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Chapman et al. (2005). These three sources are not detected at > 3.5a in the 'super-map' but 
they are detected in SCUBA photometry observations with a lower statistical significance. 

Our 13 SMG IRS targets have 850 /im and 24 /im flux densities of 2-10 mJy and 200- 
1200 /iJy, respectively. This sample is expected to be representative of the submm population 
discovered in existing blank-field 850 /im imaging. There is no evidence supporting the idea 
that SMGs which are fainter than 200 /iJy at 24 /im are fundamentally different besides lying 
at slightly higher redshifts (e.g. Pope et al. 2006). We list the positions and mid-IR fluxes for 
these 13 SMGs in Table [U SMG IDs starting with 'GN' and 'C' are from the SCUBA super- 
map (Pope et al. 2006; Pope 2007) and Chapman et al. (2005), respectively, and the 850 /im 
fluxes can be found in these papers. The majority of the SMGs are undetected in the deep 
Spitzer 70 /im images (see also Huynh et al. 2007). 11/13 SMGs have optical spectroscopic 
redshifts while the other two have photometric redshifts. The redshift distribution of our 
sample is consistent with that of all SMGs (Chapman et al. 2005; Pope et al. 2006). Fig. [1] 
shows multi- wavelength images, which span roughly the width of the IRS slit, of our targets. 
In some cases there are several optical or IRAC sources in the aperture. For GN04, GN07 and 
GN19, the two IRAC sources are both expected to be associated with the submm emission 
since both components are at the same redshift (see Pope et al. 2006). The IRS spectra of 
these 3 SMGs will have contributions from both components. For the remaining sources, 
the 8, 16 and 24 /im images show that only our primary targets are producing the mid-IR 
emission which we will detect with the IRS. 
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Fig. 1. — Images of SMG targets, in order of increasing redshift. Images are 10" x 10" 
(roughly the size of the aperture of the IRS observations) at B, i, 3.6, 8.0, 16 and 24 fim 
from left to right. Data are not available in some of these wavebands for source C3 since it 
is outside the HST GOODS-N coverage. GN04, GN07, GN19 and GN39a and b, have two 
sources within the slit and which contribute flux in the IRS spectra (See text for details). 
For the other sources, the 8, 16 and 24 /iin images show that only the primary targets are 
bright in the mid-IR. 
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Fig. 1. — Continued. 
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The Spitzer IRS observations are taken in spectral staring mode which observed the 
target at two nod positions within the slit. We observed in low resolution (R = A/AA ~100) 
using the Short-Low 1 (SL1; 7.4-14.5 /xm), Long-Low 2 (LL2; 14.0-21.3 /mi) and/or Long 
Low 1 (LL1; 19.5-38.0 jum) orders, depending on the redshift of the source. We chose the 
wavelength coverage in order to include one or more of the main PAH features at 6.2, 7.7 
and 11.3 /xm. Integration times varied depending on the flux of the target at 24 /im. We 
used blue peak-up acquisition on nearby, isolated, bright 2MASS stars. Our observations 
were taken in April and May 2006. 

Table [2] shows our list of targets and integration times. 

3. Data analysis 

We started with the S14.0.0 Spitzer IRS pipeline datcQ which produces 2-dimensional 
(2D) Basic Calibrated Data (BCD) files. Note that our targets are not resolved with Spitzer 
at these wavelengths and therefore they can be treated as point sources. 

At mid-IR wavelengths, the dominant source of noise is the sky background, mostly 
from zodiacal light. Since we have long integrations of very faint targets, additional data 
reduction steps are required to alleviate the detector and sky noise. The first step is to 
identify and clean the rogue pixels. There are a number of hot pixels in the IRS arrays 
(~ 15%). Rogue pixel masks were created from the campaign masks and from examining 
each 2D spectral image individually. We used IRSCLEANg to replace the values of the hot 
pixels by extrapolating from the surrounding pixels. Note that IRSCLEAN fails when there 
is a group of several hot pixels together, so we made sure to exclude these from the spectral 
extraction. The fraction of pixels in the 2D files which remain unusable after IRSCLEAN is 
only ~ 2%. 

The next step is to remove the latent charge build-up on the array. In long IRS in- 
tegrations, it has been found that latent charge builds up with time, despite the detectors 
being reset at the end of each integration (see Teplitz et al. 2007). This build-up is differ- 
ent depending on the wavelength (y position on the 2D spectral image We found charge 
build-up on the arrays for integrations of more than 1 hour in all LL observations. For the 
SL1 observations, we only saw the latent charge build-up for a few targets, even though 



^ttp : //ssc . spitzer . caltech . edu/ irs/ dh/ 

2 http : / / ssc . spitzer . caltech . edu/ archanaly/ contributed/irsclean /IRSCLEAN_MASK . html 
3 http : / / ssc . Spitzer . caltech . edu/ irs/documents/ irs_ultra deepjnemo .pdf 
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all of these observations were greater than 1 hour. The latent charge build-up was fit by a 
polynomial and the temporal component subtracted off. 

We explore several methods for removing the sky background including: a) subtracting 
the average spectrum at each nod position; b) subtracting individual observations from 
opposite nod positions and then averaging; c) subtracting a 'supersky' created from all files 
from that Astronomical Observation Request (AOR); and d) subtracting a supersky created 
from all our data taken in a given campaign. In all methods, we mask out all sources which 
fall within the IRS slit identified from the deep 24 /zm imaging of the GOODS-N field. The 
supersky subtraction was carried out with and without normalizing to the median in each 
image. The third method of subtracting a normalized supersky created from all observations 
from that AOR gave the lowest residual sky noise (and thus highest SNR) in the final 
spectrum therefore this is what we used to remove the sky from all our observations. 

After the sky was removed from the individual data files, we then coadded them using 
a clipped median. The coadds were checked again for persisting hot pixels and cleaned once 
more using IRSCLEAN. These steps produced two 2D spectra for each target, one for each 
nod position. 

We used the SPitzer IRS Custom Extraction (SPicrfl) to extract the ID spectra at 
each nod position. We used a narrow extraction window of 2 pixels (~10 arcseconds for LL 
observations) which was constant as a function of wavelength. We found this produced higher 
SNR spectra than using an extraction window of increasing width, starting with 2 pixels at 
the blue end. Using the optimal extraction option in SPICE, we found comparable results to 
the 2 pixel, fixed width, extraction. For each science spectrum, we also extracted a residual 
sky spectrum, offset from any sources, to provide a level of the error in our final spectrum as 
a function of wavelength. The uncertainty is calculated as the standard deviation of 10 pixels 
surrounding each wavelength in the residual sky spectrum. For each source, the ID spectra 
at each nod position were then averaged together to produce the final spectrum for each 
source. Further details about the final noise in our IRS spectra are given in Appendix |A] 

In order to calibrate the final IRS spectra, we reduced standard star calibrator observa- 
tions from the same campaign as our observations, in exactly the same way that we reduced 
the science observations. We then extracted the spectrum for the calibrator using our 2 
pixel width and the default SPICE width in order to determine the calibration factor as a 
function of wavelength. This calibration was verified against the broad-band photometry at 
8, 16 and 24 jum and we found very good agreement (see Fig. [2]). 



4 http : / / ssc . spitzer . caltech . edu/postbcd/spice .html 
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Fig. 2. — Calibration of IRS spectra compared to the MIPS and peak-up photometry. MIPS 
photometry is from the GOODS 24 /im imaging (Chary et al. in preparation) and the 16 /iin 
peak-up imaging is from Teplitz et al. (2005). Most of the 16 /im IRS fluxes are systematically 
higher because roughly half of the 16 /xm filter falls outside of the LL2 IRS coverage and so 
the IRS flux is overestimated, since the spectrum actually drops off at lower wavelengths. 
These sources are indicated as the open circles. For sources with SL1 and LL2 observations, 
there should be no such discrepancy, which we confirm. The solid point in the right panel 
which is significantly above the line is GN39b whose IRS spectrum is known to contain a 
contribution (1/3 of the flux) from GN39a. 
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3.1. Spectral line measurements 

In order to measure the individual PAH line luminosities we used the IDEA (isap|) 
environment in the SMART0 software package (Higdon et al. 2004). For each of the 4 
main PAH lines (6.2, 7.7, 8.6 and 11.3 /mi), we isolated a region around each line and 
simultaneously fit the line and continuum. Both the center and the width of the line were 
allowed to vary. All measured line widths were consistent within the errors with the line 
widths measured in IRS spectra of local star forming galaxies (Smith et al. 2007). Note that 
the errors we calculated from the residual sky were propagated in the fitting of individual 
lines. All SMG line measurements used in this paper including 3cr upper limits when the 
line is undetected in the IRS spectrum are listed in Table 

The biggest uncertainty in comparing PAH line luminosities from different samples 
comes from differences in how the spectral lines are measured, more specifically, how the 
continuum is estimated for each line. This can lead to variations in the line flux and equiv- 
alent width of up to a factor of 4 (see Sajina et al. 2007; Smith et al. 2007). Our method 
of measuring the line strengths gives similar results to the cubic spline continuum fits per- 
formed in Brandl et al. (2006, see next paragraph). This method is known to underestimate 
the flux of the PAH lines (particularly the 7.7 /mi line) compared to methods which fit all 
PAH lines simultaneously (e.g. Sajina et al. 2007; Smith et al. 2007) since the continuum 
we adopt is higher. Therefore, we took care to measure the PAH lines in any comparison 
samples in the exact same way as we did for the SMGs. 

In this paper, we compare the SMG PAH line luminosities to those of other galaxy 
populations. Brandl et al. (2006) presented IRS spectra of 22 local starburst galaxies with 
infrared luminosities ranging from 5 x 10 9 -5 x 1O 11 L . These galaxies are known to have 
little or no AGN contributing to their infrared luminosity. Because the choice for estimating 
the continuum under the PAH lines makes a significant difference to the results, we have 
measured the lines for the starburst galaxies in exactly the same way as we did for the 
SMGs. Another interesting comparison sample is the local ULIRG galaxies. Armus et 
al. (2007) presented IRS spectroscopy of 10 ULIRGs from the Bright Galaxy Sample (BGS, 
5*60 > 5.4 mJy, Soifer et al. 1987) and the IRS spectra for a sample of over 100 ULIRGs 
was recently presented in Desai et al. (2007) and Farrah et al. (2007). The local ULIRGs 
show a broad range of mid-IR properties and a mix of AGN and starburst mid-IR emission. 



5 The ISO Spectral Analysis Package (ISAP) is a joint development by the LWS and SWS Instrument 
Teams and Data Centers. Contributing institutes are CESR, IAS, IPAC, MPE, RAL and SRON. 

6 SMART was developed by the IRS Team at Cornell University and is available through the Spitzer 
Science Center at Caltech. 
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However, we did not have access to the IRS spectra for the full sample of local ULIRGs in 
order to measure the line luminosities using our method therefore we could not include them 
in this paper. 

We measured the PAH lines in 10 of the Brandl et al. (2006) starburst galaxies to verify 
that our method is consistent with their measurements. On average the luminosities of the 
6.2 and 11.3 fim lines we measure are consistent with the values listed in Brandl et al. (2006), 
but the 7.7^m line luminosity is 1.5 times larger on average from our measurements. This is 
probably because we are using a straight line for the continuum as opposed to a spline. We 
have corrected the measurements of the 7.7 /zm line luminosity for the Brandl et al. (2006) 
starbursts for this factor of 1.5 for all plots in this paper. For the Brandl et al. (2006) 
starbursts we also applied a correction to the line measurements to account for the fractional 
flux found within the IRS apertures (see Table 2 of Brandl et al. 2006) when comparing to 
the unresolved SMGs. The error on the line fluxes for the Brandl et al. (2006) starbursts is 
dominated by the calibration error which is on the order of 10% (less than the size of the 
symbols in most of the plots in this paper). 



4. Results 

In this section we use the mid-IR spectra to explore the nature of SMGs. We start by 
using the identification of PAH lines to estimate redshifts. We then proceed to decompose 
the mid-IR spectra into starburst and AGN components; starting first with a composite 
SMG spectrum since the SNR is higher and then proceeding to decompose each individual 
SMG spectrum. We then use the information from the spectral decomposition to classify the 
fractional contribution of AGN to the mid-IR luminosity and compare this with other AGN 
indicators such as optical spectroscopy and X-ray emission. These results are then used to 
guide the fitting of the full infrared SED to constrain the total infrared luminosity in the 
Section 5. 

Fig. E] shows the rest-frame (unsmoothed) IRS spectra for the 13 SMGs. PAH emission 
lines are present in all sources. Despite the relatively low SNR in the spectra, all SMGs 
seem to show very similar mid-IR spectral shapes with the exception of CI. While this 
SMG clearly shows PAH emission, it is superimposed on a steeply rising continuum. We will 
discuss this source further in Section [B] 
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Fig. 3. — IRS spectra of SMGs. The solid black curves are the raw (unsmoothed) IRS spectra, 
while the shaded region is the associated la noise from the sky background. The red dashed 
curves show the best-fit SED which is made up of an extincted power-law component (blue 
dashed, consistent with zero in most cases) and an extincted starburst component (green 
dashed). Recall that the extinction curve is not monotonic in wavelength and contains 
various features, the most notable being the 9.7/im silicate feature. 
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Fig. 3. — Continued. 




Fig. 3. — Continued. 
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4.1. Redshifts 

Given the presence of PAH features in all of our sources we are able to extract redshifts 
from the IRS spectra. In order to calculate the redshift and its associated error from the IRS 
data we follow the following steps. We assumed the following rest-frame wavelengths for the 
main PAH features; 6.22, 7.71, 8.61 and 11.27 /jm (Draine & Li 2007). For each PAH line, we 
find the center of the Gaussian fit to the line from SMART and calculate the redshift. The 
redshift for each source is calculated by averaging the redshift estimates from each PAH line. 
The total error in the redshift is the quadratic sum of the deviation due to the spread in the 
redshift from the different PAH lines and the centroid uncertainty. For sources whose spectra 
only contain 1 PAH line, the error is just the centroid uncertainty of that line. The deviation 
in the redshift from each line is calculated as la = ^/Yli( z i ~ ( z )) 2 / y/n—1. We verified 
the centroid uncertainties by performing Monte Carlo simulations. We take a Gaussian of 
the same height and width as each PAH feature and we run 1000 trials where we add noise 
which has the same standard deviation as the noise in the real data. For each trial, we fit 
the data to get the line center and keep track of the difference between the input line center 
and the measured line center. The centroid error is estimated from the standard deviation 
of this difference for all 1000 trials. We repeat this for each source. With the exception of 
GN26 whose spectrum has a very high SNR, the centroid uncertainties for individual PAH 
features range from 0.08-0.25 /xm which corresponds to 0.5-1.5 pixels. This is roughly an 
order of magnitude higher than what we would expect from the instrument based on the 
pixel size in the case of no noise. The total error in the IRS redshift ranges from 0.01-0.05 
for these SMGs. 
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Fig. 4. — Comparison of IRS spectroscopic redshifts and previous optical redshifts. All opti- 
cal redshifts are spectroscopic expect for the two sources indicated by open symbols. There 
appear to be 2 optical spectroscopic redshifts which disagree with our new IRS redshifts. In 
these cases we believe that the IRS redshift is correct for the SMG, since the mid-IR emission 
is more likely to trace the submm emission, than the optical/UV emission. 
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Table [3] lists the new IRS redshifts and the previous optical, spectroscopic or photo- 
metric, redshifts. The new IRS mid-IR redshifts are compared with the optical redshifts 
in Fig. HI Two of our 13 sources seem to have very different redshifts as determined from 
IRS and optical spectroscopy (C2 and GN06) which we discuss in Appendix [B]). Since the 
mid-IR and submm wavelengths both trace dust emission, we believe that the IRS redshifts 
presented here are correct for the galaxies which are causing the submm emission. 



4.2. Spectral decomposition of SMG composite spectrum 

We can study the ensemble properties of the SMG sample by averaging together the 
individual mid-IR spectra. When plotted in rest-frame luminosity units, the spectra for 
all SMGs in this sample show very little dynamic range. Therefore we are able to average 
the spectra without having to normalize each at a specific wavelength which can introduce 
additional biases in the composite spectrum. We exclude CI from the composite spectrum, 
since it shows a very different spectral shape from all the other SMGs (see Section |B]). The 
error on the composite spectrum is calculated from combining the errors on the individual 
flux measurements. We restrict the final SMG composite spectrum to areas where there are 
> 3 data files present at that wavelength which results in a wavelength coverage of ~5- 
12 /iin in the rest frame. Fig. shows our composite spectrum of 12 SMGs calculated from 
averaging the individual spectra (excluding CI). 
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Fig. 5. — Composite IRS spectrum of 12 SMGs (excluding CI). This SMG composite 
spectrum (black curve) has been smoothed by a Gaussian with a full width at half maximum 
(FWHM) of 0.12 /mi (rest-frame) which corresponds to the average resolution of the IRS. 
The shaded region shows the la uncertainty. Left panel: The red curve is the best-fit model 
(Equation [T]) which is composed of a PAH template (green curve) and a power- law continuum 
component (blue curve), each with extinction applied to them. Right panel: The composite 
is fit using individual PAH line profiles and a continuum, each with extinction. In either 
case, we find that the PAH emission dominates the mid-IR luminosity. 
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We detect the 6.2, 7.7 and 11.3 /im emission lines, with a hint of the 8.6 /im emission 
line in the composite spectrum. The composite spectrum also shows a very flat underlying 
continuum. In order to quantify the fractional contribution of star-formation and AGN 
activity to the emission at these wavelengths, we fit our composite spectrum with a model 
containing the 3 components discussed in Section [TJ The PAH emission is fit using two 
different templates: 1) the mid-IR spectrum of M82 (Forster Schreiber et al. 2003) and, 
2) the starburst composite template from Brandl et al. (2006). M82 is a prototypical star 
forming galaxy and its mid-IR spectrum is expected to be dominated by emission from star- 
forming regions. The starburst template from Brandl et al. (2006) is a composite of 13 local 
starburst galaxies with bright flux densities and without a strong AGN component. As we 
show below, these templates are composed almost entirely of the sum of the individual PAH 
features. The AGN emission is characterized by a power law with both the normalization 
and slope as free parameters. For the extinction, we obtain r v from the Draine (2003) 
extinction curves and applied it to both the PAH and continuum components separately. 
The extinction curve is not just monotonic in wavelength and contains silicate absorption 
features, the most notable being at 9.7 /im. Our model, F u , can be expressed as 

F v = Cl v- C2 e-° iT » + c 4 te- C5 T " (1) 

where f v is the PAH template. We performed a x 2 minimization fit for each of the two PAH 
templates and solved simultaneously for ci, c 2 , c 3 , c 4 and c 5 . 

The best fit model, F u , is shown in the left panel of Fig. [5] as the red dashed curve. The 
blue dash-dot and green dotted curves show the individual contributions from the continuum 
and PAH components, respectively. The best-fit PAH template was the starburst composite 
from Brandl et al. (2006) although M82 gave very similar results. With this PAH template, 
the extinction needed is r g . 7 ~ 1. The slope of the power law is very shallow (c 2 ~ 0.2). The 
continuum component accounts for 30% of the luminosity at these wavelengths (5-11.5 /im). 
This implies that the mid-IR luminosity of SMGs is dominated by star formation with a 
maximum mid-IR AGN contribution of 30%. 

To be clear about the origin of the emission in the PAH template, we also fit the SMG 
composite spectrum to individual PAH line profiles and a continuum with extinction. We 
include the PAH lines at 5.70, 6.22, 7.42, 7.60, 7.85, 8.33, 8.61, 11.23, 11.33 /im, each modeled 
as a Drude profile (e.g. Draine & Li 2007) where the height of each line is treated as a free 
parameter. The resulting best-fit model shown in the right panel of Fig. is very similar to 
the best-fit model which includes a PAH template and continuum (left panel). 
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Fig. 6. — SMG composite spectrum compared to those other ULIRGs. The SMG composite 
spectrum is the black curve and the shaded region shows the la error. The colored his- 
tograms are high redshift Spitzer 24 /zm-selected ULIRGs observed with the IRS from Sajina 
et al. (2007, S07); the green and red curves are composite spectra of the weak-PAH sources 
for 79.7 > 1 and r 9 . 7 < 1, respectively, and the blue curve is the composite spectrum for 
the strong-PAH sources. The light blue curve is the local ULIRG Arp220 (Charmandaris et 
al. 1999; Fischer et al. 1999). All curves have been normalized at ~ 7/im. The numbers in 
the legend in brackets indicate the number of sources in each composite spectrum. 



-23- 



In Fig. El we show the SMG IRS composite spectrum compared to the IRS spectra 
of z ~ 2 Spitzer 24 /im-selected ULIRGs from Sajina et al. (2007) and the local ULIRG 
Arp220. While Arp220 is often considered the typical local analogue to high redshift SMGs, 
it appears to have more silicate absorption and less 6.2 and 11.3 /im PAH emission than the 
average SMG. Recall that the sample in Sajina et al. (2007) consists of red 24 /zm-selected 
galaxies with S24 > 0.9 mJy, which is much brighter than most of the SMG population. 
The majority of the ULIRGs in the Spitzer-selected sample have weak PAH emission and 
are AGN-dominated; the green and red curves in Fig. [6] are composite spectra of the weak- 
PAH (i.e. AGN-dominated) Spitzer 24 /im-selected ULIRGs with T9.7 > 1 and rg.7 < 1, 
respectively. The blue curve is the composite spectrum for the strong- PAH Spitzer 24 /im- 
selected ULIRGs and is very similar to the SMG composite. In particular the relative line 
strengths are consistent. The strong- PAH Spitzer-selected sources represent only a small 
fraction of the Sajina et al. (2007) sample and these sources tend to have a higher detection 
rate in the (sub)mm than the rest of their sample (Sajina et al. in preparation). From Yan 
et al. (2007) and Sajina et al. (2007), the surface density of the z ~ 2 Spitzer 24 //m-selected 
ULIRGs is 10.4 deg -2 and the strong-PAH sub-sample of these sources has a surface density 
of 2.6 deg -2 . The surface density of SMGs with S 850 > 4mJy (i.e. L m > 5 x 1O 12 L , Pope et 
al. 2006) is 844 deg -2 (Coppin et al. 2006). While most submm-selected galaxies are much 
fainter at 24 /im, the strong-PAH Spitzer 24/im-selected ULIRGs appear to be a very small 
sub-sample of the SMG population. 



4.3. Spectral decomposition of individual galaxies 

Given the good agreement between our SMG composite spectrum and the spectrum of 
starburst galaxies, we proceed to fit the individual SMG mid-IR spectra to the model given 
in Equation [TJ For each of the 13 SMGs we solve for the five parameters in Equation [1] and 
plot the best-fit model in Fig. [3] as the red dashed curve. The separate starburst and AGN 
components including extinction are shown as the green and blue curves. For most SMGs, 
the continuum component is zero and not plotted in Fig. [3J We note that the spectrum 
of GN26 is even lower than the model at long wavelengths which confirms that the AGN 
component must be negligible in this source. 

From our model fits, we derived the fraction of the mid-IR luminosity which comes from 
the continuum component. The continuum component provides an upper limit to the AGN 
contribution since hot dust present in the most energetic HII regions can also contribute 
to the continuum. When the best-fit model does not include a continuum component, we 
calculate an upper limit by conservatively assuming that all of the emission around 6 /j,m is 
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coming from a continuum component. We used this fraction to classify the mid-IR spectrum 
and the results are listed in Table HI 8/13 SMGs are clearly starburst-dominated sources 
(continuum component < 50% of the mid-IR luminosity) and only 2/13 SMGs are AGN- 
dominated (continuum component > 50% of the mid-IR luminosity) in the mid-IR. The 
remaining sources have a combination of PAH and continuum emission in the mid-IR. The 
results for the individual galaxies are in agreement with the composite spectrum which 
showed a 30% AGN contribution on average. 

4.4. AGN classification 

Alexander et al. (2005) showed that deep X-ray observations are a powerful tool for 
identifying the presence of AGN in SMGs. They conclude that the majority of SMGs host 
an AGN, but that the bolometric output from SMGs is dominated by star formation because 
of the low ratio of X-ray luminosity to IR luminosity. They also stress that some fraction of 
SMGs will be AGN-dominated. Following the same method as Alexander et al. (2005) we 
list the X-ray classification in Table H] where a hard X-ray detection tells us that the X-ray 
emission is dominated by AGN activity. 

We find that the X-ray classification gives a higher fraction of AGN sources (6/13) than 
the mid-IR spectra (2/13) for SMGs with S 2i > 200 fj,Jy. For the X-ray classified AGN 
sources, the X-ray emission expected from star formation (e.g. Bauer et al. 2002) is less 
than that observed indicating that there is an AGN present. Our IRS observations show that 
this AGN is not important to the mid-IR luminosity. On the other hand, the line widths 
and line ratios in the UV/optical spectra more often indicate that these SMGs are starburst 
sources (Swinbank et al. 2004; Chapman et al. 2005) and appear to be insensitive to the 
AGN activity in many sources. 

CI is classified as a starburst from the X-ray data but shows a very strong AGN in the 
mid-IR. This SMG is likely to host a very obscured AGN, the kind of system that can help 
explain the unresolved portion of the hard X-ray background (e.g. Brandt & Hasinger 2005; 
Daddi et al. 2007). See Appendix IB1 for more discussion about CI. Combining the X-ray 
imaging with mid-IR spectroscopy provides the best census of AGN in SMGs. 

In Pope et al. (2006), we found that most SMGs showed evidence for a stellar bump in 
the 4 IRAC channels. We found only 4 sources which had a good fit to a simple power-law 
from IRAC+24 /mi photometry. Only one of these is in our IRS sample: GN04, since the 
others are too faint for spectroscopy at 24 /im. We find that the IRAC SED, the mid-IR 
spectrum and X-ray classification all converge on the conclusion that there is an energetically 
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important AGN in GN04. In contrast, the UV/optical spectra of this SMG did not reveal 
the presence of an AGN in it. 
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Fig. 7. — Spitzer color-color diagram as an AGN diagnostic (See also Ivison et al. 2004). 
The orange dashed and blue dotted curves show the positions of Mrk231 (an AGN) and M82 
(a starburst), respectively, as a function of redshift (redshift is indicated by the numbers 
along the track). The large red dots are the SMGs from this study for which we have IRS 
spectra. The 2 SMGs (CI and GN04) which have a > 50% contribution from continuum 
(AGN) emission to the mid-IR luminosity based on the IRS spectra (see Table H]) both lie off 
the M82 starburst sequence. Based on our IRS results, SMG which lie in the light shaded 
region are starburst (SB)-dominated in the mid-IR. The crosses indicate sources which are 
classified as X-ray AGN (Alexander et al. 2005). The smaller dots are the rest of the SMGs 
from Pope et al. (2006). Using our IRS results, the fraction of all SMGs which are starburst- 
dominated in the mid-IR (i.e. are found in the shaded region) is > 80%. The small open 
triangles are the Spitzer 24 /zm-selected ULIRGs from Sajina et al. (2007) with z > 1; most 
of these sources are outside the starburst region defined for SMGs. 



-27- 



Another approach for separating sources which contain a mid-IR AGN component is 
through the use of Spitzer color-color diagrams. Ivison et al. (2004) proposed that the Sg/ S4.5 
vs S 2 4,/Sg color plane could be used to identify AGN-dominated SMGs. The IRS spectra 
have allowed us to decompose the starburst and AGN components in SMGs therefore we 
are in a position to test the plot as a diagnostic. In Fig. [7J we plot this Spitzer color-color 
diagram. The orange dashed and blue dotted curves show the positions of Mrk231 (mid- 
IR spectrum from Rigopoulou et al. (1999) spliced with a fit to the Infrared Astronomical 
Satellite, IRAS, photometry) and M82 (Forster Schreiber et al. 2003), respectively, as a 
function of redshift (redshift is indicated by the numbers along the track). Based on these 
broad-band colors, AGN-dominated sources are expected to lie along the Mrk231 tracks while 
starburst-dominated systems should have similar colors to M82. Sajina, Lacy & Scott (2005) 
showed that while AGN-dominated sources do lie on this AGN sequence, there are also just 
as many high redshift PAH-dominated sources with high extinction in this region. The large 
red dots are the SMGs from this study for which we have IRS spectra. Interestingly, both 
SMGs (CI and GN04) which have a > 50% contribution from continuum (AGN) emission to 
the mid-IR luminosity based on the IRS spectra (see Table H]) lie off the starburst sequence. 
Based on our IRS results, SMG which lie in the light shaded region of Fig. [7] are starburst- 
dominated in the mid-IR. SMGs outside the shaded region likely have a combination of 
starburst and AGN emission in the mid-IR or lie at the highest redshifts (z > 4). 

This division between starburst and AGN-dominated sources based on IRS spectra 
agrees with the broad-band colors. The crosses indicate sources which are classified as X-ray 
AGN (Alexander et al. 2005). Many of these lie along the starburst sequence and have IRS 
spectra which are starburst-dominated. While the X-ray emission may be telling us that 
an AGN is present in the system, it does not seem to be a very good diagnostic of whether 
the AGN contributes to the bolometric luminosity. Sources which have AGN-dominated 
mid-IR spectra, like CI, clearly show very extreme colors in this diagram however there 
are not enough of these sources to characterize their allowable positions in this color-color 
space. The smaller dots are the rest of the SMGs from the GOODS-N submm super-map 
from Pope et al. (2006). Extrapolating our IRS results, the fraction of all SMGs which are 
starburst-dominated in the mid-IR is ~ 80%. The small open triangles are the z > 1 bright 
24/im ULIRGs from Sajina et al. (2007); very few of these sources lie within the starburst 
region defined for SMGs indicating that they have a wider range of mid-IR spectra. 



-28- 



5. Full SED fits 

In Pope et al. (2006), we fit the 24 jum, 850 /zm and radio photometry to a suite of 
models to determine the total infrared luminosity and dust temperatures. We had used the 
Chary & Elbaz (2001, hereafter CE01) templates and modified them by adding additional 
extinction from the Draine (2003) models. The CE01 templates were developed to be repre- 
sentative of local galaxies therefore they contain an inherent luminosity-temperature (L-T) 
relation. When fitting, we had allowed the models to scale so that we could solve for both 
the luminosity and the temperature independently (i.e. a model with high luminosity but 
cool average dust temperature is allowed). We had found that SMGs fit best to scaled-up 
versions of lower luminosity templates which had cooler average dust temperatures than 
local ULIRGs. In those fits the short wavelength part of the spectrum was constrained only 
by the 24 /xm flux. From the present study it is apparent that the 24 /im flux of SMGs will 
vary considerably with redshift, due to the presence of the PAH lines. Therefore, we have 
repeated the full SED fitting now that we have the mid-IR spectra. 
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Fig. 8. — IR SED of source C2. This figure demonstrates our template fitting to the full SED 
from mid-IR to radio wavelengths. The green dotted curve is the CE01 template which best 
matches the submm flux while the dot-dash dark blue curve is the CE01 template which best 
matches the 16 and 24 /zm fluxes. Both of these models are clearly poor fits to the other data 
points and have a higher total infrared luminosity than the best-fit to all the available data. 
The light blue dash-dot-dot curve is the best fit to all the data points > 7/im rest frame 
where we have allowed the CE01 templates to scale. While this provides a good fit to the 
mid-IR spectrum and the submm flux, the radio flux and the total infrared luminosity are 
underestimated. The red dashed curve is a scaled CE01 template with additional extinction 
from the Draine (2003) extinction curve which provides the best-fit to our data. We require 
both a scaling of local templates (in this case CE01) and additional extinction to fit the SED 
of high redshift SMGs. 
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Fig. [8] demonstrates why we need to allow the CE01 templates to scale (i.e. luminosity 
and temperature are both variables) and the need for an additional extinction term. Here 
we plot the mid-IR to radio SED for one of our SMG sources, C2. First we simply plot 
the CE01 template which best matches the sub mm flux only as the green dotted curve and 
the CE01 template which best matches the 16 and 24 /xm fluxes as the dot-dash dark blue 
curve. For these two models we emphasize that the CE01 templates have not been allowed to 
scale therefore they follow the local luminosity-temperature relation. Because SMGs are so 
luminous at mid-IR and submm wavelengths, the green and dark blue models have a warm 
average dust temperature (~ 50K). The light blue dash-dot-dot curve is the best fit to all 
the photometry between rest-frame 7 — 300 /xm where we have allowed the CE01 templates 
to scale. While this provides a good fit to the mid-IR spectrum and the submm flux, the 
radio flux and Ltr are underestimated compared to our best fit. We know that the CE01 
templates need to be scaled up to match the photometry of SMGs however as we scale them 
up we are introducing a greater dust mass which must then cause more extinction. The red 
dashed curve is a scaled CE01 template with additional extinction from the Draine (2003) 
extinction curve which provides the best-fit to our data. We require both a scaling of local 
templates (in this case CE01) and additional extinction to fit the SED of high redshift SMGs. 

We thus fit the IRS spectrum, the 70 /im, submm and radio photometry to the CE01 
templates plus Draine (2003) extinction for each of the 13 SMGs in this papeiu. We solve for 
the three free parameters (luminosity, average dust temperature and extinction) and find the 
template with the lowest x 2 - Fig. shows the full SED fits. In general, the resulting fits are 
very good. The dark curve and dots are the points used in the fit, while the open symbols 
show the other photometry points (mostly 16 and 24 /xm) which were not used in the fits 
but agree well with the best-fit SED. We required a scaling of 10-100 in luminosity to fit the 
local CE01 templates to the SMG photometry. Since high redshift SMGs can be well fit by 
modified CE01 templates which contain an inherent radio-IR correlation, this suggests that 
SMGs also follow a similar radio-IR correlation to local galaxies (see also Kovacs et al. 2006). 

We integrated the full SEDs in Fig. [9]from 8-1000 /xm to get the total infrared luminosity. 
Values of Ljr are listed in Tableland they are in good agreement with the values from Pope 
et al. (2006). 



7 For the GN39 double system, we split the total submm flux according to the radio flux of each component 
(note that this is the same as splitting by 24 /xm flux) . 
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Fig. 9. — Mid-IR to radio SEDs of SMGs. The IRS spectrum (dark curve) plus 70 /xm, 
submm and radio photometry (solid circles) are used in the fit. The open squares indicate 
the photometry at 16 and 24 fim. The red dashed line is the best fit modified CE01 template 
while the blue dash-dot curve is Mrk231 normalized to the AGN fraction of the mid-IR given 
in Table HI We used the latter to determine the contribution of the AGN component to Lir. 
The horizontal dotted lines indicate the region within which Lir is calculated. 
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Fig. 9. — Continued. 



27 



26 



LO 
10 

Q 25 
Q 24 
Q 23 

o 22 
o 27 

LO 26 

25 



10 

LO 

LO 

LO 
LO 

LO 

LO 

LO 

LO 

LO 
LO 

LO 

LO 25 

LO 24 

LO 23 

LO 22 
LO 27 

LO 26 



24 



23 



26 



25 



24 



23 



26 



-4 


1 1 \ 

1 1 \ 


GN06 " 




-H 1 1 — | 1 1 — 

*s 


1 — 1 — 1 1 — 1 — 1 h 

CI 


— . ' 


-H 1 1 — 1 1 1 — 

*». 

w \ 

*\ 

\ 


1 1 | 1 1 | h 

GN05 

— i — 1 — i — 1- 




-H 1 1 — 1 1 1 — 

/ 9 % 


GN19 " 

s _ — — " 




-H 1 1 — I 1 1 — 

^ N 
~/ S, V 


GN04 - 
• 

*"*•■ — .— *"■■ """ "" 



25 



24 



23 



22 



10 1 



10 2 10 3 10 4 

Rest wavelength (jim) 



10 5 



Fig. 9. 



Continued. 



-34- 



By simply scaling to the wavelength where the SED is a maximum as an estimate of 
the average dust temperature, we obtain a median of 32 K for these 13 SMGs. This is 
consistent with previous estimates of dust temperatures in SMGs (Chapman et al. 2005; 
Kovacs et al. 2006; Pope et al. 2006). We also found a similar trend between luminosity 
and temperature as found in Chapman et al. (2005) which is offset to cooler temperatures 
from the local correlation. Knowing the shape of the far-IR SED and the total luminosity, 
we can also estimate the dust mass - following the relation given in Dunne et al. (2000), the 
median dust mass for this sample is 2.4 x 1O 8 M . This is roughly an order of magnitude 
higher than dust masses observed in galaxies locally although the local galaxies studied are 
generally much less luminous. 

In order to quantify the effect of the extinction term on the models, we calculate Z^r 
without the extinction. The median amount that the extinction term changes Lir by is 13% 
(0-20%), and so models which do not include the extinction term will incorrectly determine 
Lir by around 13%. We find an average value of rg.7 ~ 1 and we don't see any correlations 
between the extinction and Lir for the SMGs given the error bars. 

For sources which we found to contain an AGN component in the mid-IR, we calculate 
the contribution to Lir from the AGN. Using ISO spectra, Tran et al. (2001) found that the 
ratio of L 5 _io Mm /LiR was about 2.5 times greater for AGN than for starbursts i.e. as you 
extrapolate from mid-IR to far-IR wavelengths AGN become less important bolometrically 
than starbursts. In order to quantify this for our sample, we scale the Mrk231 (an AGN 
whose Lir is known to be dominated by AGN emission and not star formation, Armus 
et al. 2007) template to the AGN fraction (listed in Table H]) of the mid-IR luminosity 
and then integrate the SED from 8-1000 /xm to get the contribution to Lir from the AGN 
component. For sources which have an upper limit to the AGN contribution, we plot the 
scaled Mrk231 in Fig. Oat this limit. The values of Lf^ N from integrating under the curves 
and = Lir — Lf^ N are given in Table [5J As seen in Fig. [HI only CI has enough AGN 
contribution in the mid-IR to make a significant contribution to Lir. Therefore the intense 
IR luminosities seen in SMGs are powered almost entirely by star formation. In Section 
we discuss how this luminosity can be converted to give us the star formation rate (SFR). 



6. PAH luminosities 

PAH line profiles have broad wings which contribute a significant fraction of the line 
flux. PAH line fluxes therefore, need to be consistently measured in a uniform way before 
comparisons between local sources and high redshift sources can be performed. We have 
measured line fluxes of all sources consistently, including the comparison sample of low 
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redshift starbursts (Brandl et al. 2006), using the technique described in Section |3~T1 

The PAH line luminosities and equivalent widths of the SMGs including 3a upper limits 
when the line is undetected are listed in Table The 6.2 /zm line is the cleanest to measure 
and use for diagnostics since it is the least affected by silicate extinction and line blending. 
We can see that Lpah,7.7 for the SMGs varies by over an order of magnitude, while Lpah,6.2 
and £pah,ii.3 have much smaller dynamic ranges of a factor of a few. While this might 
seem odd, we note that only a subset (40-60%) of our sample have 6.2 and 11.3 /xm spectral 
coverage whereas 92% have coverage of the 7.7/zm emission line. Furthermore, the redshift 
range for the sources with 6.2 and 11.3 fim spectral coverage is very small and therefore we 
expect a small dynamic range in PAH luminosity. The sources which are dominated by an 
AGN in the mid-IR, as classified by our IRS spectra (CI and GN04), both show suppressed 
6.2 /im emission relative to their 7.7/zm emission. 

Fig. [TUl shows the relative strengths of PAH lines compared to local starbursts (Brandl et 
al. 2006). The top panels show Lpah,7.7 as a function of Lpah,6.2 (left) and Z/pah,ii.3 (right). 
The smaller bottom panels show the residual scatter after removing the fit shown in the 
upper panels; the two larger colored error bars quantify the la scatter seen in the starbursts 
and SMGs. The outlying source in the left panel of Fig. [TO] is CI, the SMG with the largest 
AGN contribution. 
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Fig. 10. — Plot comparing the PAH line luminosities of SMGs presented in this paper (red 
circles) to that of local starbursts (blue squares, Brandl et al. 2006). All line measurements 
have been made in a consistent way (see Section I3TTT) . la error bars are shows on all points 
(the 10% errors on the starbursts are similar in size to the symbols). The top panels show 
£pah,7.7 as a function of Z/pah,6.2 (left) and Lpah,ii.3 (right) where the dashed lines shows 
the best fit to the starburst galaxies only. The similar relationship seen in starbursts and 
SMGs between PAH line strengths, indicates that SMGs are starburst-dominated and that 
PAH line strengths do not vary significantly with redshift. The smaller bottom panels show 
the residual scatter after removing the fit shown in the upper panels; the larger colored error 
bars quantify the la scatter seen in the starbursts and SMGs, respectively. 
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While the SMGs do not show enough dynamic range in PAH luminosity to pick out any 
strong correlations within the population, they do fall on the relation between various PAH 
lines established for local starburst galaxies. The dashed lines in Fig. [TD] are the best-fits to 
the local starbursts only (blue squares). We find essentially the same relation if we fit the 
local starbursts and the SMGs, indicating that the PAH luminosities in SMGs are simply 
scaled up (albeit by several orders of magnitude) from less luminous starburst galaxies. The 
best-fit relations from fitting both the SMGs and the local starburst galaxies are 



The SMG CI appears to lie off this relation with an excess of 7.7 /im PAH emission 
for a given Lpah,6.2 or Lpah,u.3- Rigopoulou et al. (1999) found suppressed 6.2 /im emission 
relative to 7.7 //m emission in ISO spectra of local ULIRGs compared to starbursts and they 
conclude that this is primarily because of extinction. This suggests that, on average, the 
SMG CI considered here may contain greater extinction or a more dominant AGN component 
than most SMGs and starburst galaxies since an AGN continuum will dilute the strength of 
the weaker 6.2 and 11.3 /zm features relative to the strong 7.7 /im line. 



log(LpAH,7. 7 ) = (1-0 ± 0.1) + (0.94 ± 0.01) l0g(i:pAH,6.2) 
log(LpAH,7.7) = (0.5 ± 0.1) + (1.00 ± 0.01) log(LpAH,11.3). 



(2) 
(3) 
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Fig. 11. — Rest-frame equivalent widths of main PAH features. Again the SMGs are denoted 
by the red circles and the local starbursts (Brandl et al. 2006) by the blue squares. All 
measurements have been made in a consistent way (see Section IBTTj) . Measurements of both 
the 6.2 and 7.7 /xm PAH lines are not available for all SMGs. CI stands out again as having 
a very small 6.2 /im EW, below the values of EWpah,6.2 = 0.2 (vertical dashed line) which is 
often used to separate AGN-dominated and starburst-dominated systems. 
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Fig.fTTl shows the equivalent widths of the 6.2 and 7.7 /zm PAH lines. The vertical dashed 
line at £Wpah,6.2 = 0.2 is often used to separate starburst-dominated and AGN-dominated 
sources (Armus et al. 2007; Sajina et al. 2007). With the exception of CI, all of our sources 
with measurable 6.2 fim PAH emission have EWpaii,6.2 > 0.2, which indicates that they 
are starburst-dominated. We have again plotted the local starbursts for comparison. We 
found that the values in Brandl et al. (2006) for EWpah,7.7 were on average 1.6 times lower 
than our measurements of the starburst EWpah,7.7> due to differences in the way the lines 
were measured. Therefore we have revised all the Brandl et al. (2006) points upwards in 
this plot. The starbursts cover a small area of this plot, while the SMGs are spread over a 
larger region, although with much larger error bars. Our current data do not have sufficient 
SNR to determine if there is any genuine spread or correlation among the EWs, however 
our error bars are sufficiently small to exclude equivalent widths of less than 0.2 in either 
PAH line (with the exception of CI). Our SMGs show similar equivalent widths to the 
PAH-dominated ULIRGs from the Spitzer 24/im selected sample in Sajina et al. (2007), if 
we decrease their values of EWparj.7 by a factor of 4 to account for the different methods 
for measuring the EWs. With -EWpah,6.2 values of > 0.5 and t 9 7 ~ 1, high redshift SMGs lie 
in region 2C of the diagnostic plot in Spoon et al. (2007). Sources in this region are classified 
as PAH-dominated and, not surprisingly, M82 also falls in this region. 

Utilizing ISO spectra of starbursts, ULIRGs and AGN, Genzel et al. (1998) used the 
ratio of the 7.7 /zm PAH line flux to the continuum flux, l/c, as a diagnostic for AGN 
contribution in the mid-IR emission. In general, galaxies with l/c > 1 are classified as 
starbursts and those with l/c < 1 are AGN-dominated. As we stress in Section [3.11 this 
type of analysis is strongly dependent on the specific choice made for the continuum fitting 
procedure. Nevertheless, we calculate l/c ~ 2.3 for the composite SMG and we note that all 
SMGs in our sample satisfy the l/c > 1 criterion for starburst-dominated sources. This is 
true even for CI which clearly shows a strongly rising continuum indicative of the presence 
of an AGN. Based on Fig. [TTJ the 6.2 /im PAH line may be more useful in weeding out 
AGN-dominated sources. 

In Fig. [12] we plot the PAH line luminosities as a function of Lip for SMGs and local 
starbursts (large panels). For the SMGs, we plot the coming from star formation 
(i.e. after subtracting the component from the AGN). Errors in Lpi for the SMGs have been 
estimated from the error in the submm flux only, since this dominates (assuming fixed CE01 
templates). In all panels, the dashed line is the best-fit to the local starburst galaxies; the 
smaller panels show the residual scatter around this best-fit relation. The best-fit relation 
for the starburst galaxies only is consistent with the fit to both the starburst galaxies and 
the SMGs. The best-fit relation for both the SMGs and starburst galaxies is given by 
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log(L PAH ,6.2) = (-2.7 ± 0.1) + (1.01 ± 0.01) log(Z IR ) 
log(LpAH,7. 7 ) = (-1.5 ± 0.1) + (0.95 ± 0.01) log(L IR ) 
log(LpAH,n.3) = (-1.6 ± 0.2) + (0.91 ± 0.01) log(L IR ). 



(4) 
(5) 
(6) 



The derived parameters are all within la of those found by fitting only the starbursts. 
The slopes are close to one indicating a direct proportionality between Lj^ and the PAH 
line luminosities. 

The first thing to note is that the SMGs lie on the relation established for local starburst 
galaxies extrapolated to very high luminosities. This is what we expect if both are dominated 
by the same emission mechanism, namely star formation. Brandl et al. (2006) showed that 
the PAH line luminosities are correlated with the IR luminosities within their starburst 
sample. Schweitzer et al. (2006) showed a similar correlation exists in a sample of QSOs 
due to star- formation in them. Peeters et al. (2004) showed a correlation between L FIR 
and I/pah,6.2 for a range of scales of star formation in local massive star-forming regions to 
normal and starburst galaxies. The fact that the SMGs lie on this relation demonstrates the 
reliability of our full SED fitting to derive the total infrared luminosity and confirms that 
SMGs do have incredibly high infrared luminosities powered primarily by star formation. 
Note that CI docs not stick out as an outlier in these plots, since we have removed the AGN 
component of its L IR . 
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Fig. 12. — Correlations between Lir and PAH luminosities. The dashed lines show the 
best-fit relations for the local starburst galaxies only; the smaller panels show the residual 
scatter around the best-fit relations. The SMGs clearly lie on the relations established by 
local starburst galaxies. 
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The second thing which stands out in Fig. [12] is that there is significant scatter in this 
relation for both the low redshift starbursts and the high redshift SMGs. In the smaller 
panels the two larger colored error bars quantify the la scatter seen in the starbursts and 
SMGs. The scatter is similar for the SMGs and starbursts and is considerably higher than 
would be implied by the measurement uncertainty. Other properties of individual objects, 
such geometry of the starburst region, spatial distribution of the dust, viewing angle relative 
to a merging event, dumpiness of the star formation, etc. are likely playing a role in this 
scatter. 

We found no correlations between the IR luminosity and equivalent widths of the PAH 
lines in SMGs. This is consistent with what was found in Brandl et al. (2006) for the local 
starburst galaxies, and shows that both the mid-IR continuum and the PAH strength scale 
together with increasing Ljr. We also find no correlation between X-ray luminosity and Lir 
(or the PAH luminosities) for the SMGs, a correlation which would have been expected had 
the IR emission been dominated by an AGN. 

Using the Kennicutt (1998) relation between L IR and SFR, we plot the corresponding 
SFR on the top axis of Fig. [T2J The Kennicutt (1998) relation between Lir and SFR assumes 
a Salpeter (1955) initial mass function and applies to starbursts with ages less than 100 Myr. 
Assuming that both the PAH and far-IR emission are tracing the star formation, we can 
use the Kennicutt (1998) relation and the above PAH line luminosity to Lir correlations to 
determine the relation between SFR and the PAH line luminosity. We use the 6.2 /im PAH 
line since it is the least affected by the silicate absorption and line blending. We find that 
given the PAH luminosity, one can calculate the SFR using: 

SFR [Moyr- 1 ] ~ lO" 7 L PAH , 6 . 2 [L ]. (7) 

Similar (although perhaps slightly less accurate) relations can be found for Lpahj.7 and 
£pah,ii.3 using Eq. H 

7. Discussion 

We have shown from the mid-IR spectra of high redshift SMGs that these systems are 
starburst-dominated and contain only a small contribution from the AGN. Despite, the fact 
that X-ray observations suggest that most SMGs harbor an AGN, only a small fraction 
of SMGs (2/13) contain an AGN which contributes significantly to the mid-IR luminosity. 
When extrapolated to the total IR luminosity, all SMGs are dominated (> 50% of L IR ) by 
star formation and not AGN activity. 
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From Table[5]we see that the median for our sample is 4x 10 12 L Q , which corresponds 
to a SFR of around 700 M yr _1 . This is an extremely high rate of star formation, capable of 
creating a massive 1O 11 M galaxy in only 10 8 years. Coupled with the low number density 
of SMGs compared to other massive galaxies at high redshift, this leads us to suspect that 
SMGs are a short-lived phase in the evolution of massive galaxies. 

SMGs are only one of several populations of high redshift galaxies which are ultralu- 
minous in the infrared. Many z ~ 2 galaxies selected in deep 24 /zm surveys also qualify as 
ULIRGs (e.g. Daddi et al. 2005; Yan et al. 2005, 2007). These other ULIRGs are generally 
not detected at submm wavelengths. What then are the fundamental differences between 
various populations of high redshift ULIRGs? 

One of the main factors which could cause a bifurcation in the population of ULIRGs at 
high redshift is the presence or absence of a bolometrically significant AGN. The evolutionary 
scenario proposed by Sanders et al. (1988) shows massive galaxies going through several main 
stages en route to becoming a massive elliptical. It starts with an IR luminous phase most 
likely triggered by a massive merger. During this stage of intense star formation, the AGN 
is also growing. As the AGN becomes larger it begins to feedback on the galaxy, eventually 
quenching the star formation completely by blowing off all the dust and gas. Thus begins the 
QSO phase, where the AGN is free to dominate the emission. Eventually the QSO consumes 
its gas and we end up with a quiescent massive elliptical galaxy. In this scenario, the submm 
emission will be at a maximum during the initial IR luminous phase (Springel, Di Matteo & 
Hernquist 2005). However, one uncertainty in this scenario is the timescale. In particular, 
how long is the IR luminous phase and how long does it take the AGN to develop? 

Since most ULIRGs at z ~ 2 show a mix of AGN and starburst activity (e.g. Houck et 
al. 2005; Sajina et al. 2007) they fit into this evolutionary scenario at the stage where the 
AGN has developed but before it has had a chance to expel all the gas and dust, since the 
IR luminosities are high. Since SMGs show less of an AGN contribution, it is likely that 
they represent a slightly earlier phase of this scenario, before the AGN has time to develop 
fully and is significant to the IR luminosity. This is also consistent with the finding that 
SMGs have cool dust temperatures (Pope et al. 2006), partly because the star-formation is 
taking place on extended spatial scales and partly because the AGN is not yet energetically 
important for heating the dust. Other high redshift ULIRGs, like those selected at 24 fim 
(e.g. Sajina et al. 2007), should have warmer dust temperatures (or have overestimated 
infrared luminosities), otherwise they would be detected in the submm. Therefore, they 
likely represent the stage after the SMG phase in this evolutionary scenario. 

For local ULIRGs, there is evidence which shows that the luminosity of CO emission 
decreases as the merger progresses (Rigopoulou et al. 1999), which is consistent with de- 
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creasing submm emission in the later stages of the merger. This idea is also consistent with 
morphological studies which show that SMGs are often tidally disrupted systems (Chapman 
et al. 2003; Conselice et al. 2003; Pope et al. 2005). Alternatively, SMGs and other high 
redshift ULIRGs might appear different because of different merger progenitors. 

Since the number density of SMGs is much less than that of all 24 /im selected IR- 
luminous galaxies, then this phase in the evolutionary sequence must be shorter. We can 
make a rough estimate for the duration of the submm luminous phase. In the 2 Gyrs between 
1.5 < z < 3 there are roughly 20 SMGs and 500 24 /zm-selected galaxies (with L IR > 1O U L , 
see Chary et al. 2004) within a portion of the GOODS-N field. Assuming that all IR-luminous 
24 /im-selected galaxies will undergo a submm luminous stage during this period of time, we 
calculate the duration of the submm luminous phase to be on the order of 10 8 yrg^|. As 
mentioned above, this timescale is what is needed to create a galaxy with a stellar mass of 
~ 1O 11 M (Borys et al. 2005) at a rate of ~ 1000M Q yr _1 . This timescale is consistent with 
those derived from gas masses of SMGs (Greve et al. 2005). 



8. Conclusions 

Spitzer IRS spectroscopy has been obtained for a sample of 13 SMGs brighter than 
200 /iJy at 24 /im. This has effectively confirmed the identification of these SMG counter- 
parts, and in some cases has provided a new redshift estimate. The SMGs show strong PAH 
emission and, on average, only a small continuum contribution from hot dust. 

We have explored several diagnostics from the mid-IR spectra of SMGs to determine 
the level of AGN contributing to the luminosity at these wavelengths. All of them seem to 
converge on a picture in which SMGs are starburst-dominated systems with at most a 30% 
contribution from an AGN at mid-IR wavelengths. Their mid-IR spectra are similar to a 
scaled spectrum of local starburst galaxies like M82. 

The classification of an SMG as AGN or starburst-dominated from the mid-IR often 
disagrees with the classification from the X-rays. While X-ray observations are better suited 
for determining the presence of an AGN in an SMG, the mid-IR spectra can determine how 
important the AGN is to the total infrared luminosity since it is directly detecting the hot 
dust emission and not subject to details about the geometry of the obscuring material. 



8 It is not clear what fraction of 24 /im-selected galaxies will undergo a submm-luminous phase since SMGs 
are known to be massive and could represent only the most extreme galaxies. Nevertheless, this calculation 
provides a useful reality check on the expected duration of the SMG phase. 
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Full IR SED fits to the mid-IR data, far-infrared and radio photometry show that SMGs 
are best-fit by scaled up versions of local IR-luminous galaxy templates with additional ex- 
tinction from the Draine (2003) extinction curves. These models have cool dust temperatures 
(T = 32 K) and high Lirs which imply SFRs of ~ 700 M yr _1 assuming a Salpeter initial 
mass function. 

SMGs lie on the local relation between L IR and Lpah,6.2 (or Lpahj.7 or Lpah,ii.3)? which 
means that the PAH line flux can be used to estimate the SFRs in these systems, albeit with 
large uncertainties. Equivalent widths, which have been the focus of some other studies, are 
much less useful in this regard. 

SMGs appear to be an earlier phase in the evolution of massive galaxies than other 
high redshift ULIRGs. The average high redshift ULIRG shows warmer dust and less PAH 
emission. It is therefore likely that the AGN in SMGs has not yet grown strong enough to 
heat the dust substantially and destroy the PAH molecules. 

We have been able to put constraints on the contribution from AGN activity to the rest- 
frame 5-12 /xm mid-IR emission, and have extrapolated this to far-IR wavelengths. Without 
more data from 30-100 /im rest frame it is difficult to directly probe the emission at far-IR 
wavelengths. We will be able to model this part of the spectrum in more detail once data 
from the Herschel Space Observator^ and SCUBA-2 (Holland et al. 2006) become available. 
Details of the SMG phase will be understood only after the spatially resolved spectroscopy 
capabilities of ALMA0 become available. 
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quality of this paper. We thank Anna Sajina for helpful discussions and for providing 
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A. Spitzer IRS noise 

Fig. [T3] shows the final RMS of our spectra as a function of the integration time and 
spectral order. Values from the Spitzer Science Center sensitivity calculator; SPEC-PET0 
are given for comparison as the dashed curves. We confirm that the sensitivity goes as 
cr(mJy) = a x t(s)~ 1//2 even for very long integrations. The scatter in the measured values is 
probably dominated by the accuracy of the sky subtraction in these crowded fields. Given the 
scatter, our observed sensitivities are consistent with those from SPEC-PET. The decrease 
in sensitivity in the SL1 order could be related to our choice of 240 second ramps in fields 
with backgrounds > 20MJysr _1 at 24 //m, which results in incorrect droop correction (see 
the Spitzer Observer's Manual). 



http : //ssc . Spitzer . caltech . edu/tools/ specpet/ 
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Fig. 13. — Observed la RMS as a function of integration time for our IRS staring mode 
observations. All observations from our Spitzer GO2-20456 program are shown here. The 
dashed curves show the predicted values from SPEC-PET for each order. We confirm that 
the sensitivity goes as a oc t -1 / 2 even for very long integrations. The scatter in the measure 
values is probably dominated by the accuracy of the sky subtraction in these crowded fields. 
Given the scatter, our observed sensitivities are consistent with those from SPEC-PET. 
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B. Notes on individual sources 

Here we describe several of our sources in more detail. Sources which are not discussed 
here are relatively straightforward in terms of interpreting their mid-IR spectra. 

B.l. GN06 

The optical and IRS redshifts disagree for this SMG. The IRS redshift is based on 
3 PAH features: 6.2, 7.7 and 11.3/im and the optical spectroscopic redshift of 1.865 is 
inconsistent with the positions of these emission features, even within the errors. The optical 
spectroscopic redshift of GN06 comes from Chapman et al. (2005) where they note that the 
optical galaxy was offset from the radio counterpart by about an arcsecond. From the 
images of GN06 in Fig. [TJ we see there are two optical knots both of which are offset from 
the counterpart which is detected in the radio and mid-IR. For this source, we conclude that 
the optical and mid-IR spectra are picking up different galaxies and the submm emission is 
coming from the mid-IR detected galaxy. 

B.2. C2 

C2 has the largest discrepancy between the IRS and UV/optical spectroscopic redshifts. 
The IRS redshift for this source is based on the 7.7 and 11.3 /im PAH features and there 
is also a hint of the 8.6 /xm line. The optical spectroscopic redshift would have placed the 
7.7 //m PAH line right in the middle of this spectrum and therefore it is inconsistent with 
the IRS spectrum. The optical spectroscopic redshift comes from detections of both Ha and 
Lymana emission (Swinbank et al. 2004; Chapman et al. 2005). The images in Fig. [T]show 
an extended optical galaxy. This is a curious source, however since the mid-IR is more likely 
to trace the submm emission, we conclude that the IRS redshift is correct for this SMG. 

B.3. CI 

This SMG is not detected in the deep Chandra 2 Msec imaging of GOODS-N. Further- 
more, the optical spectrum from Swinbank et al. (2004) clearly shows Ha, Nil and the SII 
doublet and the low Nil to Ha ratio is used to classify this SMG as a starburst. However 
the mid-IR spectrum of CI clearly shows a steeply rising continuum indicative of an AGN 
emission. This is also apparent from the extreme Spitzer colors shown in Fig. [71 The mid-IR 
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spectrum also shows PAH emission and the ratio of PAH luminosity to IR luminosity for 
CI is significantly lower than those seen in local starbursts. We conclude that this SMG 
harbors an Compton-thick AGN which accounts for 44% of the total IR luminosity. This 
is the only SMG which contains a bolometrically significant AGN. Interestingly, CI has not 
been detected in CO (Greve et al. 2005). From the X-ray flux limit, we can place a limit on 
the column density of Ah ^> 10 24 cm~ 2 . CI is relatively bright at 24 //m (1.2 mJy) suggesting 
that these types of obscured AGN at high redshift should be detectable in even shallow 
Spitzer surveys. 



B.4. GN39 double system 

One further SMG deserves special discussion and that is GN39. This SMG has two ra- 
dio/MIPS counterparts separated by 8 arcseconds, which we refer to as GN39a and GN39b. 
Both are confirmed to lie at z — 1.996 (Swinbank et al. 2004; Chapman et al. 2005), there- 
fore it is likely that they are both contributing to the submm emission. The separation 
corresponds to 70 kpc at this redshift. 

Our IRS observations picked up both counterparts and, with a narrow extraction win- 
dow, we were able to extract a separate spectrum of each component. Some of the flux 
from GN39a falls into the spectral extraction window of GN39b however a comparison of 
the 24 fxm flux from the IRS spectra to that from the MIPS 24 fim image confirms that 
only 1/3 of the flux in the IRS spectrum of GN39b is from GN39a. The extracted spectra 
confirm that these two sources are indeed at the same redshift and both spectra look very 
similar to each other and to the rest of the SMGs, supporting the idea that they are both 
submm emitters. Furthermore, both components are X-ray detected and classified as ob- 
scured AGN by Alexander et al. (2005), although the two IRS spectra show a completely 
starburst-dominated system at mid-IR wavelengths. 

Since these two SMGs are at the same redshift, we know that they must be associated. 
Given the physical distance between them, larger than those seen in most local ULIRGs 
(Murphy et al. 1996), we could be seeing the early stage in the merger of two massive 
galaxies (Springel, Di Matteo & Hernquist 2005). Note that GN04, GN07 and GN19 also 
have double radio counterparts, which are thought to be at the same redshift. However, in 
these sources, the two components are separated by only 3 arcseconds and therefore they do 
not resolve into two separate 24 /im sources. Thus the IRS spectra for these sources contain 
contributions from both components. 
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Table 1. Coordinates and flux densities of IRS targets 



SMG Name 


SMG ID 


RA 




DEC 






Sl6 


•524 




SVo 


SMMJ.. 




J2000 




J2000 






(mjy) 


(mjy) 


(mjy) 


123653.1+621120 


GN31 


12:36:53. 


22 


62: 


11:16 


.7 





.301 


± 0.007 


0.367 


± 


0.006 


< 3.5 


123635.5+621238 


GN26 


12:36:34. 


51 


62: 


12:40 


.9 


0. 


.992 


± 0.012 


0.446 


± 


0.005 


13.9 ± 1.8 


123701.2+621147 


GN17 


12:37:01. 


59 


62: 


11:46 


.2 


0. 


.209 


± 0.006 


0.710 


± 


0.008 


< 3.5 


123622.6+621629 


C2 


12:36:22. 


66 


62: 


16:29 


.5 





.098 


± 0.006 


0.414 


± 


0.007 


< 3.5 


123555.1+620901 


C3 


12:35:55. 


13 


62: 


09:01. 


.6 






n/a 


0.374 


± 


0.009 


< 3.5 


123711.1+621325 


GN39a 


12:37:11. 


37 


62: 


13:31 


.1 


0. 


.097 


± 0.005 


0.537 


± 


0.009 


< 3.5 


123711.9+621331 


GN39b 


12:37:11. 


97 


62: 


13:25 


.8 





.051 


± 0.005 


0.225 


± 


0.007 


< 3.5 


123621.3+621711 


GN07 


12:36:21. 


27 


62: 


17:08 


.1 





.048 


± 0.007 


0.370 


± 


0.011 


< 3.5 


123618.7+621553 


GN06 


12:36:18. 


33 


62: 


15:50 


.4 


0. 


036 


± 0.006 


0.330 


± 


0.008 


< 3.5 


123600.2+621047 


CI 


12:36:00. 


16 


62: 


10:47 


.3 


0. 


.478 


± 0.011 


1.220 


± 


0.008 


5.5 ± 1.2 


123618.8+621008 


GN05 


12:36:19. 


13 


62: 


10:04 


.3 


0. 


.044 


± 0.007 


0.215 


± 


0.006 


< 3.5 


123707.7+621411 


GN19 


12:37:07. 


19 


62: 


14:08 


.0 






< 0.02 


0.280 


± 


0.015 


< 3.5 


123616.6+621520 


GN04 


12:36:16. 


11 


62: 


15:13 


.5 





060 


± 0.006 


0.303 


± 


0.007 


< 3.5 
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Table 2. Spitzer IRS observations. 



SMG ID 


SLl(x240s) 


Integration time 
LL2(xl20s) LLl(xl20s) 


Total (hr) 


GN31 




35x2 




2.3 


GN26 


15x2 


20x2 


[12x2] 


3.3 [4.1] 


GN17 


20x2 


12x2 




3.5 


C2 






45x2 


3.0 


C3 




31x2 


40x2 


4.7 


GN39a 




20x2 


28x2 


3.2 


GN39b 




[20x2] a 


[28x2] 


[3.2] 


GN07 




35x2 


45x2 


5.3 


GN06 




35x2 


55x2 


6.0 


CI 




6x2 


10x2 


1.1 


GN05 






55x2 


3.7 


GN19 




[20x2] 


[28x2] 


[3.2] 


GN04 






70x2 


4.7 



a Sources with time in square brackets are secondary sources observed in the slits of the 
primary targets. 
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Table 3. Measured IRS redshifts compared to optical spectroscopic and photometric 

redshifts. 



SMG ID 


% mag (AB) a 


~ d 
^optical 


Z mS ± 1(7 


Comment 


GN31 


21.8 


0.935 h 


0.93 ± 0.03 




GN26 


22.7 


1.219 6 


1.23 ± 0.01 




GN17 


27.7 


1.72^ 


1.73 ± 0.01 


new IRS redshift 


C2 


25.5 


2.466 e ' / 


1.79 ± 0.04 


inconsistent 


C3 


23.5 6 


1.875 6 


1.88 ± 0.02 




GN39a 


25.9 


1.996' 


1.98 ± 0.01 




GN39b 


26.4 


1.992 6 


1.99 ± 0.04 




GN07 


27.8/23.9 c 


1.988 6 , 1.992' 


1.99 ± 0.02 




GN06 


27.8 


1.865 6 


2.00 ± 0.03 


inconsistent 


CI 


24.8 


1.994 6 , 2.002' 


2.01 ± 0.05 




GN05 


24.9 


2.60 9 


2.21 ± 0.03 


new IRS redshift 


GN19 


25.4/>28 c 


2.484 e , 2.490' 


2.48 ± 0.03 




GN04 


26.2/>28 c 


2.578 e 


2.55 ± 0.01 





a i magnitudes are from the GOODS ACS observations unless otherwise noted. 
b Capak et al. (2004). 

c These SMGs have two counterparts within the IRS aperture. 

d Optical redshifts are spectroscopic unless otherwise noted. 

c Chapman et al. (2005). 

f Swinbank et al. (2004). 

g Pope et al. (2006), photometric redshift. 

h Cowie et al. (2004). 
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Table 4. Classification of mid-IR spectra compared with classifications based on X-ray 

imaging and UV/optical spectroscopy. 



SMG ID 




Classification" 


Mid-IR 




X-ray ft 


UV/Op c 


Mid-IR 


continuum% d 


GN31 


U 


n/a 


SB 


< 45 


GN26 


SB 


SB 


SB 


< 15 


GN17 


SB 


n/a 


SB 


< 29 


C2 


AGN 


SB 


SB 


< 34 


C3 


AGN 


SB 


SB+AGN 


48 


GN39a 


oAGN 


SB 


SB 


< 35 


GN39b 


oAGN 


SB 


SB 


< 10 


GN07 


SB 


SB 


SB+wcak AGN 


18 


GN06 


SB 


SB 


SB+AGN 


47 


CI 


SB 


SB 


AGN+wcak SB 


82 


GN05 


U 


n/a 


SB 


< 7 


GN19 


oAGN 


SB 


SB 


< 14 


GN04 


AGN 


SB 


AGN+SB 


61 



a oAGN = obscured AGN, SB=starburst, U=undetected. 

b Classified as AGN if detected in the hard X-ray band. 

c Classification from Chapman et al. (2005) and/or Swinbank et al. (2004). 

d Perccntage of mid-IR flux that comes from continuum emission. 
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Table 5. PAH strengths and infrared luminosities. 



SMG ID 


PAH luminosity (1O 9 L ) 


PAH equivalent width' 


' (Mm) 


L IR (10 12 L ( 






6.2 fim 


7.7 


11.3 fim 


6.2 fim 


7.7 fim 


11.3 fim 


SB 


AGN 


Total 


GN31 




3.4 ± 1.2 






0.32 ±0.11 




> 0.17 


< 0.10 


0.27 


GN26 


12.2 ± 1.2 


42.9 ±5.1 


8.1 ± 1.4 


0.38 ±0.04 


1.00 ±0.13 


1.18 ±0.54 


> 3.5 


< 0.1 


3.6 


GN17 


13.0 ±4.1 






0.22 ±0.07 






> 1.3 


< 0.7 


2.0 


C2 




31.1 ±6.3 


7.0 ±1.8 




1.08 ±0.30 


1.00 ±0.29 


> 2.8 


< 0.5 


3.3 


C3 


< 19.2 


9.2 ±2.7 


7.9 ±2.7 


< 0.50 


0.25 ±0.08 


0.45 ±0.16 


5.2 


0.8 


6.0 


GN39a 


13.9 ±4.9 


39.2 ± 12.9 


< 36.3 


0.66 ±0.26 


0.86 ±0.30 


< 0.86 


> 4.2 


< 0.8 


5.0 


GN39b 


12.3 ± 5.6 


23.1 ±5.1 


< 29.4 


1.09 ±0.65 


0.59 ±0.14 


< 1.14 


> 2.7 


< 0.1 


2.8 


GN07 


12.5 ± 2.6 


17.6 ±3.2 


9.9 ± 1.8 


1.10 ±0.38 


0.41 ±0.08 


1.58 ±0.47 


6.7 


0.2 


6.9 


GN06 


18.0 ±1.5 


26.3 ±4.4 


10.5 ± 2.2 


0.83 ±0.10 


0.93 ±0.16 


2.95 ± 1.38 


5.7 


0.7 


6.4 


CI 


6.3 ± 1.5 


59.1 ± 12.6 


< 18.4 


0.05 ±0.01 


0.54 ±0.14 


< 0.18 


5.4 


4.2 


9.6 


GN05 




12.9 ±4.3 






0.44 ±0.15 




> 3.9 


< 0.1 


4.0 


GN19 


12.5 ± 3.0 


32.5 ±9.0 




0.59 ±0.20 


0.77 ±0.22 




> 11.9 


< 0.4 


12.3 


GN04 


< 23.7 


43.3 ±9.9 




< 0.45 


0.63 ±0.15 




4.5 


2.1 


6.6 



a The equivalent widths are all in rest frame. 



